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nonlinear. In nonlinear, diffusive shock acceleration, compression ratios will be higher and the shocked temperature 
lower than test-particle, Rankine-Hugoniot relations predict. Furthermore, the heating of the gas to X-ray emitting 
temperatures is strongly coupled to the acceleration of cosmic-ray electrons and ions, thus nonlinear processes which 
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on X-ray line models. Here we apply an algebraic model of nonlinear acceleration, combined with SNR evolution, to 
model the radio and X-ray continuum of Kepler's SNR. 
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INTRODUCTION 

More than twenty years of spacecraft observations in 
the heliosphere have proven that collisionless shocks can 
accelerate particles with high efficiency, i.e., 10-50% of 
the ram energy can go into superthermal particles (e.g., 
Eichler [17]; Gosling et al. [21]; Ellison et al. [19]; Tera- 
sawa et al. [30]). Energetic particles exist throughout the 
universe and shocks are commonly associated with them, 
confirming that shock acceleration is important beyond 
the heliosphere as well. In fact, shocks in supernova rem- 
nants (SNRs) are believed to be the main source of Galac- 
tic cosmic rays, and these shocks are expected to be much 
stronger than those in the heliosphere and can only be 
more efficient and nonlinear. 

The conjecture that collisionless shocks are efficient 
accelerators is strengthened by results from plasma sim- 
ulations, which show efficient shock acceleration consis- 
tent with spacecraft observations (e.g., Scholer, Trattner, 
& Kucharek [29]; Giacalone et al. [20]), and other indi- 
rect evidence comes from radio emission from SNRs (see 
Reynolds & Ellison [27]) and equipartition arguments in 
AGNs and y-ray bursts (see Blandford & Eichler [6] for 
an early review). There is also clear evidence that shocks 
can produce strong self-generated turbulence. This has 
long been seen in heliospheric shocks (e.g., Lee [24, 25]; 
Kennel et al. [23]; Baring et al. [3]) and there is evidence 
that it occurs at SNRs as well (i.e., Achterberg, Bland- 
ford, & Reynolds [1]). 



While the importance of nonlinear (NL) shock acceler- 
ation is evident, NL solutions to diffusive shock acceler- 
ation are complicated and results are often unwieldy and 
difficult to use for astrophysical applications. Therefore, 
we have developed a simple, algebraic model of diffu- 
sive shock acceleration, based on more complete studies, 
which includes the essential nonlinear effects (Berezhko 
& Ellison [4]; Ellison, Berezhko, & Baring [18]). This 
technique is computationally fast and easy-to-use, yet 
includes (i) the modification of particle spectra when 
the backpressure from energetic ions smooths the shock 
structure, and (ii) the influences on the shock dynamics 
when the magnetic turbulence is strongly amplified by 
wave-particle interactions. 

The complications of NL shock acceleration and the 
many parameters required to characterize it are offset 
somewhat by the fact that the entire particle distribution 
function, from thermal to the highest energies, is inter- 
connected and must be accounted for self-consistently 
with a nonthermal tail connecting the quasi-thermal pop- 
ulation to the energetic one. Because energy is conserved, 
a change in the production efficiency of the highest en- 
ergy particles must impact the thermal properties of the 
shock heated gas and vice versa. If more energy goes into 
relativistic particles, less is available to heat the gas. In 
contrast, the power laws assumed by test-particle models 
have no connection with the thermal gas, energy conser- 
vation does not constrain the normalization of the power 
law, and the spectral index can be changed with no feed- 
back on the thermal plasma. Furthermore, there is a direct 



linkage between protons and electrons (which produce 
most of the photon emission associated with shocked gas) 
in nonlinear models, so the entire emission from radio to 
gamma-rays, plus cosmic -ray observations, can, in prin- 
ciple, be used to constrain the models. 

Here we describe some of the nonlinear features ex- 
pected to occur in young SNRs and investigate some 
implications of efficient cosmic -ray production on the 
broad-band continuum from Kepler's SNR. We refer to 
Berezhko & Ellison [4] and Ellison, Berezhko, & Baring 
[18] for details of the NL shock model and its applica- 
tion to particle and photon production in SNRs. Work on 
the NL X-ray line emission from Kepler is in progress, 
i.e., Decourchelle, Ellison, & Ballet [12]. Previous test- 
particle calculations of the X-ray emission from Kepler 
have been reported by Borkowski, Sarazin, & Blondin 
[7], who used a two-dimensional hydrodynamic simu- 
lation, and by Rothenfiug et al. [28], who investigated 
the emission from the reverse shock. The work of De- 
courchelle & Ballet [11] is mentioned below. 

NONLINEAR SHOCK ACCELERATION 

The nonlinear effects in shock acceleration are of two 
basic kinds: (i) the self-generation of magnetic turbulence 
by accelerated particles and (ii) the modification (i.e., 
smoothing) of the shock structure by the backpressure of 
accelerated particles. Briefly, (i) occurs when counter- 
streaming accelerated particles produce turbulence in the 
upstream magnetic field which amplifies as it is convected 
through the shock. This amplified turbulence results in 
stronger scattering of the particles, and hence to more ac- 
celeration, quickly leading to saturated turbulence levels 
near 8B/B ~ 1 in strong shocks. The wave-particle in- 
teractions produce heating in the shock precursor which 
may be observable. Effect (ii) results in the overall com- 
pression ratio, r to t, being an ever increasing function of 
Mach number (as if the effective ratio of specific heats 
Yeff — * 1 an d analogous to radiative shocks), i.e., 

r tot ~ 1.3M^ 4 if Ml >M A0 , (1) 

or by 

r tot ~1.5M^ 8 (2) 

in the opposite case (Mso is the sonic and Mao is the 
Alfven Mach number). Simultaneously, shock smooth- 
ing causes the viscous subshock to be weak (r su b <C r tot ) 
and the temperature of the shocked gas to drop below the 
test-particle value. 

There are many parameters associated with NL shock 
acceleration (see Ellison, Berezhko, & Baring [18] for a 
listing), but the most important ones that determine the 
solution are the Mach numbers (i.e., the shock speed, 
pre-shock density, and magnetic field) and the injection 
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FIGURE 1. Comparison of nonlinear and test-particle shock 
results as a function of upstream sonic Mach number. The top 
two panels show that, at Mso > 100, the downstream density 
can be > 5 times larger and the shocked temperature > 10 times 
smaller if the acceleration is NL. The bottom panel shows that 
for Mso > 100 and T\ m j p = 10~ 3 , more than 70% of the pressure 
in the shocked gas ends up in relativistic particles in NL shock 
acceleration. These results depend on parameters other than 
Mso, an d in particular are a strong function of r|j n j p , particularly 
at high Mso- The heavy-weight curves are for r|; n j p = 10~ 3 and 
the light-weight curves are for r|; n j p = 10 -4 . 



efficiency, T| i n j p (i.e., the fraction of total protons which 
end up with superthermal energies). As described in 
Berezhko & Ellison ([4]), we use Alfven heating in the 
precursor which reduces the efficiency compared to adia- 
batic heating. Significantly, parameters typical of young 
SNRs should result in NL acceleration. 

Figure 1 shows a comparison between test-particle 
(dashed lines) and nonlinear results as a function of sonic 
Mach number, Mso- As illustrated in the top panel, the 
overall compression ratio, r to t, is generally > 4 and can 



increase without limit in the NL case for r|i n j p = 10 . 
Simultaneously, the subshock compression ratio, r su b, is 
less than 4 (typically r su b ~ 3). The overall compression 
determines the shocked density, while the subshock is 
mainly responsible for heating the gas. Thus, the shocked 
gas will be cooler (middle panel) and denser if the accel- 
eration is NL compared to a test-particle (TP) case. This 
is an important consideration for X-ray line models. The 
particle spectrum at the highest energies is determined by 
r tot , and thus will be flatter than in a TP shock. How- 
ever, the spectral shape at lower energies depends on r su b, 
indicating that the particle spectrum will go from being 
steeper than the TP prediction at low energies to flatter at 
the highest energies, this is the so-called concave signa- 
ture of NL shock acceleration. The bottom panel shows 
that a large fraction of the shocked pressure can end up in 
relativistic particles (dotted line) for large Mso- These re- 
sults depend on the injection efficiency assumed and the 
turnover in r tot at large Mso seen in the r|i n j,p = 10~ 4 ex- 
ample is a transition to a strong, unmodified shock (see 
Berezhko & Ellison [4] for a full discussion). 

In Figure 2 we show the phase space distribution func- 
tions, f(p), for electrons (dashed line), protons (solid 
line), and helium (dotted line). The parameters used for 
this forward shock are those discussed below for Ke- 
pler's SNR and given in Table 1. These particles pro- 
duce photons from synchrotron, bremsstrahlung, inverse- 
Compton, and pion-decay processes (see Baring et al. 
[2] for details) and these spectra are shown in Figure 3 
for both the forward and reverse shocks. As discussed 
below, we have chosen parameters to cause the reverse 
shock to contribute substantially to the total X-ray emis- 
sion, but this choice impacts the emission over the en- 
tire range from radio to gamma-rays. While the radio 
synchrotron comes primarily from the reverse shock, the 
forward shock dominates emission at TeV energies with 
a pion-decay component. Any changes in parameters to 
accommodate the X-ray observations must be consistent 
with observational constraints in all other bands. 

APPLICATION TO KEPLER'S SNR 

We now apply our model to the radio and X-ray emis- 
sion seen in Kepler's SNR. The X-ray emission shows the 
presence of lines and these have been interpreted, in mod- 
els without acceleration (see Decourchelle & Petre [13]; 
Decourchelle et al.. [14]), as coming from the ejecta ma- 
terial heated by the reverse shock. However, these pre- 
vious models also required the presence of a power law 
photon continuum, presumably coming from synchrotron 
emitting TeV electrons (e.g., Reynolds [26]). 

Here, we use the evolutionary model of Truelove & 
McKee [31] to obtain the forward and reverse shock pa- 
rameters at the ~ 395 yr age of Kepler's SNR and, us- 
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FIGURE 2. Downstream phase space distribution functions, /, 
versus momentum, p. We have multiplied f(p) by [p/(m p c)] 4 
to flatten the spectra, and by [(m p c) 3 /n p o] to make them di- 
mensionless (ii p o is the upstream proton number density). The 
solid line shows protons, the dotted line shows helium, and the 
dashed line shows electrons. 
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FIGURE 3. Photon spectra from the forward (FS) and reverse 
shocks (RS) obtained with parameters used to generate the spec- 
tra shown in Figure 2. Notice that emission in the X-ray band 
(1-100 keV) is totally bremsstrahlung for the FS but has a syn- 
chrotron component in the RS. Emission at gamma-ray energies 
is dominated by the FS, while the radio comes totally from the 
RS. 



ing these shock parameters, calculate the NL particle and 
then photon spectra. Figure 4 shows the synchrotron and 
bremsstrahlung emission from Figure 3 and indicates that 
we are able to obtain a good fit to the data within the 
constraint that the continuum emission from the reverse 
shock contribute substantially to the X-ray emission (our 
model only produces continuum which lies below any X- 
ray emission lines). Since the forward and reverse shocks 
have very different emission profiles (Figure 3), the abil- 
ity to resolve the individual shocks in radio and X-rays 
will be exceptionally important for constraining the mod- 
els. 

We note that the relation between the line and con- 
tinuum emission is not straightforward and requires a 
detailed calculation as is currently in progress, i.e., De- 
courchelle, Ellison, & Ballet [12]. This is particularly the 
case in young SNRs because, while the lines and con- 
tinuum both depend on the electronic temperature in the 
shocked gas, the lines also depend strongly on the ejecta 
composition and on the ionization state. The ionization 
state can be strongly influenced by non-equilibrium ef- 
fects and by the history of the shocked plasma. What can 
be said, however, is that the presence of lines restricts the 
synchrotron intensity to some level comparable to (with a 
factor of 10 say) or below the bremsstrahlung continuum. 

As mentioned above, NL shocks are extremely com- 
plicated with many (often poorly defined) parameters and 
our purpose in this preliminary work is only to give some 
indication of the issues involved when NL effects are con- 
sidered. The basic point is that, if nonlinear cosmic-ray 
production occurs, the X-ray modeling cannot be done 
without considering emission in other bands, particularly 
radio and gamma-rays (if available). For Kepler we know 
that strong lines are present, indicating that the thermal 
gas behind the reverse shock is contributing substantially 
to the emission. This strongly constrains the range of ac- 
ceptable parameters. Previous TP models have also indi- 
cated that a continuum component is likely to be present. 
If so, this is likely to be synchrotron from TeV electrons 
and must be consistent with the radio observations. 

The parameters used to generate Figure 4 are given in 
Table 1 . All terms in this Table are explained in Ellison, 
Berezhko, & Baring ([18]) but some additional explana- 
tion is required since we include here, for the first time, a 
calculation of acceleration at the reverse shock, (i) The 
unshocked proton number density, n p o, is arbitrarily cho- 
sen for the FS and determined for the RS with eqs. (20), 
(28), and (30) in Truelove & McKee [31]. (ii) The un- 
shocked magnetic field, Bo, is arbitrarily chosen and the 
shocked field (which produces the synchrotron emission) 
is taken to be B2 = r tot Bo- We allow for an independent 
Bo for the reverse shock and, in fact, the only way the 
RS can contribute substantially to the X-ray emission is 
for the RS Bq to be much larger (25 vs. 1.5 /jG) than 
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FIGURE 4. Radio and X-ray data for Kepler's SNR compared 
to the emission from the forward and reverse shocks. The radio 
data is from Reynolds & Ellison [27] and X-ray data is adapted 
from Decourchelle & Petre [13] (Decourchelle, private commu- 
nication). We use a standard vF v representation and have sup- 
pressed the error bars on the X-ray points. 



the FS Bq. (iii) The supernova energy, E sn , and ejecta 
mass, M e j, are standard values, (iv) We assume the same 
unshocked proton temperature, T p o = 10 4 K, for the FS 
and RS. The solutions are relatively insensitive to this pa- 
rameter, (v) The shocked electron to proton temperature 
ratio, r e 2/r p 2, is a free parameter in the NL shock model 
and is set to 1 for both shocks, (vi) The injection effi- 
ciency, r|i n j,p, is an important model parameter determin- 
ing the overall acceleration efficiency. Values much larger 
than rjinj p = 2 x 10~ 4 tend to produce strong nonthermal 
tails on the bremsstrahlung emission which may be incon- 
sistent with the observations. Values much smaller than 
2x 10 -4 yield test-particle solutions with typical temper- 
atures higher than allowed by the X-ray observations. We 
note that adjustments in T|j n j p produce extremely large 
modifications in the TeV results and can be highly con- 
strained with gamma-ray observations, (vii) The elec- 
tron to proton ratio at relativistic energies, (e/p) le \, is an 
arbitrary parameter in the simple model but is expected 
to be between 0.01 and 0.05 if (e/p) re \ in galactic cos- 
mic rays is typical of that produced by the strong shocks 
in young SNRs. This factor is very important for deter- 
mining the pion-decay contribution to TeV gamma-rays. 
Furthermore, lowering (e/ p) re \ lowers the overall emis- 
sion in the radio to X-ray band but increases the relative 
contribution of thermal to synchrotron photons in the X- 
ray band. A firm determination of the synchrotron contri- 
bution to X-rays, combined with a gamma-ray detection, 
will help determine (e/p) re i, one of the most important 



unknown parameters in shock acceleration, (viii) The 
maximum energy cosmic rays obtain in the NL model 
depends on the scattering mean free path, X, which is as- 
sumed to be, 

"X = T| m fp '"g.max (fg/ fg,max) , (3) 

where r| m f p is taken to be independent of particle momen- 
tum (Baring et al. [2], use the notation T|), r s = p/(qB) is 
the gyroradius in SI units, r g . max is the gyroradius at the 
maximum momentum, p max , and a is a constant param- 
eter (a = 1 = T| m fp is roughly the Bohm limit). Here we 
only consider a = 1. A fairly large value of r| m fp = 15 
(giving a maximum cosmic -ray energy well below that 
obtained in the Bohm limit) is required to avoid hav- 
ing the X-ray emission totally dominated by synchrotron 
from TeV electrons, (ix) We assume that the initial den- 
sity profile in the ejecta has a power law density distribu- 
tion, p oc r~", with n = 9, and that the unshocked ISM is 
uniform. 

The output values are also given in Table 1 . Of these, 
Vsk and R s k are taken directly from the Truelove & McKee 
[31]) solution. Once these are obtained, the Mach num- 
bers are determined (given the relevant input parameters), 
and then the NL shock model determines the compres- 
sion ratios, temperatures, and particle spectra. Using the 
particle spectra, the continuum emission is calculated. It 
is important to notice that the shocks are highly efficient 
and NL. They yield total compression ratios > 4, shocked 
temperatures nearly 10 times lower than the correspond- 
ing test-particle shocks, and place the majority of the to- 
tal energy flux, e re i, into relativistic particles (mainly pro- 
tons), yielding y e ff's less than 5/3. 

Finally, our model includes a rough estimate of the 
emission volume 

Venus « (4z/3)R 3 Jr m , (4) 

as described in Ellison, Berezhko, & Baring [18]. This 
volume is considerably less than the total remnant vol- 
ume. If we assume a distance to Kepler of D sm = 5 
kpc (e.g., Decourchelle & Ballet [11]), the normalization 
given by eq. (4), and shown in the Table, matches the ob- 
servations. 

DISCUSSION AND CONCLUSIONS 

X-ray line and continuum emission contains a vast 
amount of information on supernova (SN) ejecta elemen- 
tal composition, ISM density and elemental composition, 
the SN explosion energy, and the mass of ejecta. In 
addition to heating the plasma, the forward and reverse 
shocks accelerate some fraction of the shocked material 
to cosmic-ray energies and this acceleration is believed 
to be quite efficient, removing energy from the thermal 



Table 1. Parameters for Kepler's SNR model 



Input parameters 


Forward shock 


Reverse shock 


r — 3n 

n p0 [cm J ] 


0.5 


1.25 




1.5 


25 


E sn [10 51 erg] 


1 




M ei \m. : 


5 




T p0 [K] 


10 4 


10 4 


?e2/7p2 


1 


1 


Tlinj.p 


2xl0~ 4 


2xl0~ 4 


(e/p)td 


0.03 


0.03 


Tlmfp 


15 


15 


/I 


9 


9 


Output values 


V sk [kms- 1 ] 


4100 


1750 


Ask [PC] 


2.5 


2.1 




260 


110 


M A o 


1040 


43 


''lot 


11 


6.4 


r.sub 


3.9 


3.5 


B 2 [jjG] 


16 


160 


^max.p 


4.1x10" 


1.7xl0 12 


^|inj,e ^ 


5.1xl0~ 4 


6.3 xl0~ 4 


T p2 [K] 


2.8 xlO 7 


1.3 xlO 7 


Tt P [K]** 


2.0 xlO 8 


9.4 xlO 7 


£ rel 


0.79 


0.54 


Yeff 


1.40 


1.45 


Flux parameters 


Am- [kpc] 


5 




Vemis [PC 3 ] 


6.0 


6.4 



* This is the maximum proton energy, but since synchrotron 
losses are unimportant here, electrons have the same £ max . 
' This is the electron injection efficiency. 

** This is the temperature the shock gas would have obtained if 
no acceleration took place. 



plasma (e.g., Kang & Jones [22]; Dorfi & Bohringer [16]; 
Berezhko, Ksenofontov, & Petukhov [5]). Despite the ex- 
pected efficiency of shock acceleration, virtually all cur- 
rent X-ray line models assume that the shocks that heat 
the gas do not place a significant fraction of their en- 
ergy in cosmic rays (exceptions to this are Chevalier [10] 
and Dorfi [15]). Here, we investigated the broad-band 
continuum emission expected in Kepler's SNR from effi- 
cient shock acceleration, by coupling self-similar hydro- 
dynamics (Chevalier [10]; Truelove & McKee [31]) with 
nonlinear diffusive shock acceleration (e.g., Berezhko, 
Ksenofontov, & Petukhov [5]). We were able to show 
that the radio and X-ray continuum can be fit with reason- 
able parameters in a way that allows the reverse shock to 



contribute substantially to the total X-ray emission. This 
constraint is required since X-ray line models of Kepler 
(e.g., Decourchelle & Petre [13]) require emission from 
the shock-heated, metal-rich ejecta material. 

This preliminary calculation is not fully self-consistent 
for several reasons. Most importantly, we use self-similar 
results (i.e., Truelove & McKee [31]) to model the SNR 
evolution. These solutions neglect the effects of energetic 
particle escape from the FS and assume that the ratio of 
cosmic -ray pressure to total pressure at the shock front is 
a constant. They also assume that y e ff = 5/3. As indicated 
in the bottom panel of Figure 1 , the total shocked pressure 
doesn't differ much with or without cosmic -ray produc- 
tion and we have demonstrated (i.e., Decourchelle, Elli- 
son, & Ballet [12]) that NL shock results do not change 
dramatically over most of the age of Kepler for typical 
values of the injection parameter, r|i n j p > 10~ 4 . How- 
ever for lower r|i n j p , the nonlinear solutions can have test- 
particle, unmodified solutions at very high sonic Mach 
numbers with a rapid transition to the NL solution as the 
Mach number drops (see Fig. 1 and Berezhko & Ellison 
[4] for a detailed discussion). The self-similar solutions 
we use are still approximate, however, because we have 
not yet modified them for the change in y e ff that results 
when a substantial fraction of the shocked pressure is in 
cosmic rays. We also neglect cosmic-ray diffusion and 
assume they are spatially coupled to the gas - an excel- 
lent approximation for all but the highest energy particles. 
However, the highest energy electrons produce the X-ray 
synchrotron photons so there may be differences that are 
not modeled in the emission volumes and other important 
parameters between the radio and X-ray bands. Finally, 
we have not included absorption in our models which is 
probably not important for the Kepler radio observations, 
but will be required to model the low energy X-rays. 1 

Besides providing a more self-consistent model of 
photon production, predictions from NL shock models 
provide a test of the fundamental assumption that SNRs 
are the primary source of galactic cosmic-ray ions. If 
this is so, the acceleration is almost certainly nonlinear 
since 5-30% of the total ejecta kinetic energy is required 
to replenish cosmic rays as they escape from the galaxy. 
Since shocks put more energy into accelerated ions than 
electrons, nonlinear effects seen in X-ray emission will 
be evidence for the efficient shock acceleration of ions as 
well as electrons. X-ray observations potentially provide 
in situ information on cosmic-ray ion production, com- 
plementing observations of pion-decay y-rays in this re- 
gard. Any inference of nonthermal tails on electron dis- 



1 Note that we have only plotted the X-ray observations above 2 keV in 
Figure 4 to avoid conflict with the absorbed low energy end of the X-ray 
distribution. 



tributions in X-ray observations will provide information 
on electron injection, the least well understood aspect 
of shock acceleration. Such in situ information on high 
Mach number shocks is available nowhere else. 

Our modeling of Kepler's SNR suggests that typical 
source parameters produce large nonlinear effects in the 
broad-band spectrum and suggest that the test-particle ap- 
proximations that are almost universally used are inad- 
equate for SNRs as young as this. Besides the differ- 
ences discussed above, one might expect that the growth 
rate of the Rayleigh-Taylor instability will be greater in a 
cosmic-ray modified shock because of the larger spatial 
gradients of density, pressure, etc. Furthermore, the high 
compression ratios result in a considerably thinner region 
between the forward and reverse shocks than predicted in 
the TP case (Decourchelle, Ellison, & Ballet [12]). This 
places the contact discontinuity closer to the shock and 
may make it easier for the Rayleigh-Taylor "fingers" to 
distort or overtake the FS, a situation that appears diffi- 
cult with normal TP parameters (Chevalier, Blondin, & 
Emmering [9]; Chevalier & Blondin [8]). Another im- 
portant difference concerns electron heating and equili- 
bration. The higher densities in the NL models mean 
that electron heating is much more efficient than in TP 
shocks. Our initial calculations [12] suggest that it may 
be possible to obtain full equipartition between electrons 
and ions, at least for high values of r|i n j p , in the shocked 
ejecta. 
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